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We derived the low energy effective action for the collective modes in asym- 
metric fermionic systems with attractive interaction. We obtained the phase 
diagram in terms of the chemical potentials. It features a stable gapless su- 
perfluidity with one Fermi surface on the BEG side of the resonance. Also we 
predict a sharp increase in outer core of a vortex, i.e. vortex size, upon entering 
into the gapless phase. This may serve as a signature of a gapless phase. 

Keywords: BCS-BEC crossover; asymmetric ultracold Fermi systems; gapless 
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1. Introduction 

Experiments in cold atomic gasses renewed the interest in fermionic pair- 
ing. In particular, much effort was made to understand superfluidity 
in imbalanced fermionic gasses. Imbalanced systems have unequal num- 
ber of particles of two fermionic species that pair. The system consists of 
fermions of two different species, which correspond to two hyperfine states 
of a fermionic atom such as ^Li. These fermions have opposite spin and 
interaction between them can be tuned by using a Feshbach resonance.'^ 

For zero imbalance, the system properties are qualitatively well under- 
stood using mean field. ^ In weak coupling the system lives in a weakly cou- 
pled BCS state and crosses over to a strongly coupled BEC state through 
the resonance region. While the extreme BCS and BEC regimes are in good 
control in mean field theory, quantitative understanding of the phases close 
to resonance comes mainly from Monte-Carlo calculations.^ At resonance, 
fluctuations change mean field only quantitatively. Situation is completely 
different at nonzero polarization, where fluctuations change mean field re- 
sults around the resonance qualitatively; and many features of the phase 
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diagram arc not caught by the mean field. For imbalanced systems, a quah- 
tatively complete picture of the phase diagram has not been established yet. 
At small asymmetry, and weak coupling, there is standard BCS superflu- 
idity, which breaks down at Clogston limit. At large asymmetry, proposed 
possibilities are phase-separation,^' gapless (or breached) superfluidity, 
deformed Fermi sea pairing, and non- homogeneous or LOFF pairing. ""^^ 

The phase diagram of the system at T=0 as a function of scattering 
length and the chemical potential difference has been explored. In the mean 
field approximation,^'^ it was found that on the BCS side of the resonance 
there are no stable gapless superfluid phases with net polarization; and 
at the resonance the first order phase transition from superfluid to nor- 
mal phase takes place as asymmetry is increased, without any intervening 
gapless superfluids. In experiment with trapped atoms this leads that the 
system phase separates, because the gapped phase can not feature a net po- 
larization, unpolarized superfluid is in the central region of the trap and a 
polarized normal fluid is at the exterior. There was an effort done to go 
beyond the mean field approximation by using Monte-Carlo simulations,^^ 
an expansion in e = d — 4 space dimensions,^'' 1/A'' loop expansion around 
the BCS-BEC solution, ^^-^^ a superfluid local density approximation. 

The authors of^^ propose a phase diagram in the plane of scaled po- 
larization T] 6n and diluteness parameter k = —1/na^ which features 
a splitting point near the resonance at nonzero mismatch where the ho- 
mogeneous superfluid, a LOFF like inhomogeneous phase, and the gapless 
superfluid phase coexist. They also find stable gapless fermionic modes with 
one and two Fermi surfaces on the BCS side of the resonance. This picture 
is supported by Monte-Carlo simulations.-'^^ 

In order to understand how fluctuations in the condensate around the 
mean field value affect the phase diagram of cold atomic gasses with asym- 
metry, we study the effective theory describing these fiuctuations. 

2. Effective low energy theory 

We consider a non-relativistic system consisting of two species of fermions 
■^1 and '^2 of equal mass m but different chemical potentials /xi and /K2- We 
assume that the Feshbach interaction between fermions of different species, 
given via the scattering length, can be modeled by a point like four Fermi 
interaction with coupling g = l/(fci?as), n = fc|./(27r^). BCS is at g ^ — oo, 
Qg < 0; BEC is at 5 — > oo, > 0. The effect of the attractive interaction 
between fermions is to produce a difermion condensate. Condensate sponta- 
neously breaks U{1) associated with conservation of total fermion number 
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(fermions arc not charged); hence there is one Goldstone mode and the 
system is superfluid. By gauging U{1), condensate spontaneously breaks 
local U{1); hence the gauge boson becomes massive, i.e. there is Meissner 
mass and the system is a superconductor. We established a correspondence 
between cocfBcicnts of effective theory for Goldstone mode and screening 
masses of the gauge field. Mean field Lagrangian is quadratic in fermion 
fields, 

^^t-^{p) + ^^ia^ -A{x)e \. |A(a; 



V A*(a;)e idt + ^{p) - Sfxa^ J A ^' 

^' stands for the four component Nambu-Gorkov spinor. We consider a 
homogeneous condensate, independent of x, i.e. A(a;) = A, as a mean-field 

solution. The quasiparticle dispersion law is e± = ±(5/z + ^'^(p)^ + A^. At 
some momenta it features zero energies, describing spherical Fermi surfaces 
in momentum space. 

We include fluctuations of the condensate around the mean field sohi- 
tion, by introducing the field rj{x), A{x) = A + r]{x). Partition function 
contains fermion field ^ and fluctuation of the condensate fields, 

z = J r),y*r)7yr)*tx)^-e-'5[*^'*''''''*i (2) 

To get an effective action for the fiuctuations, we integrate out the fermion 
fields; we get the fermion determinator, 

s[v,v*] = j d^x[^\A + n{x)\''] 

1 2 -(A + r;*(,T)) -^^.^^{p)+5^i)i 

We expand the fermion determinator in fluctuations and their derivatives, 
S[r],r]*] = S^^^ -\-S^^^ +<S^^^ + ... . As a result we get an effective theory of 
collective modes (Goldstone and Higgs) away from unitarity. To the leading 
order. S^'^'^ is the free energy of the system without fluctuations. Due to the 
gap equation the next term 5^^-' vanishes (stationary point of the action). 
The lowest nontrivial term in the expansion is 5^^) , 

c(2)__l?^VrAf fcW ,../Mn(fc)Mi2(fc)WA(fc)\ 

with rj = {\ + iO)/ \pl. It is UV finite; UV divergent contributions cancel 

exactly by employing the gap equation. (Note, the gap equation is UV 
divergent and it is regulated by the cutoff A, usually the Fermi momentum. 
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which is the highest scale in the problem). We separate magnitude and the 
phase of the fluctuations, 77 = (A + i0)/\/2. The physical meaning of the 
real and complex components of the fluctuation field is easy to understand 
in small fluctuation and long wavelength limit, 



/^^[A{dtcl>{x)f - -{dicl>{x)f] (6) 



-C + Dkl-^k'^ ikoF 
-ikoF Ak^ - f 

Let us first consider the phase of condensate, but not its magnitude, 
fluctuates, A — > Ae"^^^^. The field represents the Goldstone mode associ- 
ated with the spontaneous breaking of the total fermion number, 711+712. 
There is no explicit breaking, hence the mass of Goldstone boson is zero. 
In the low energy Lagrangian density, 

B. 
3 

A is the kinetic energy, B is the spatial variation of the Goldstone mode. 

Negative A OT B means instability towards the growth of phase fluctuations. 
If A and B are positive, the system is stable and the speed of sound is 
^/B/3A. In our analyses we flnd ^ > always, it can be B < 0. 

The case when magnitude fluctuates corresponds to the Higgs mode, 
A ^ A + X{x)/y/2, with the Lagrangian of a massive bosonic field 

= -IcXixf + lDidtX{x)f - ^{d,X{x)f (7) 

with all positive coefficients, C, D, E > 0. = C/D is the mass squared, 
squared of the gap in the excitation spectrum. CA^ is change in free energy 
5^°) due to change in magnitude of condensate; C is proportional to the 
curvature of the potential, C > is local minimum (stable/meta-stable) 
and C < is local maximum (unstable). When D < or E < 0, mean field 
solution is unstable with respect to time or space-dependent fluctuations of 
the magnitude; D > always. 

B < 0, E < mean growth of spatially non-uniform fluctuations, unsta- 
ble to phase and magnitude modulation of the condensate. However, these 
instabilities point to different ground state of the condensate. B < 0: con- 
densate develops phase modulation which carries a current, balanced by a 
counter-propagating current carried by gapless fermions. E < 0: formation 
of a spatial modulation in the magnitude of the condensate, which does not 
carry a current. 

F mixes Higgs and Goldstone components. F vanishes in weak coupling. 
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Fig. 1. Left panel. The {5fi, fi) plane is divided into three regions corresponding to 
different number of gapless surfaces in momentum space, marked 0, I, II (solid lines). 
Unstable regions according to criteria 1 and 3: speed of sound for Goldstone and Higgs 
modes should be real (local criteria), free energy of the superfluid phase is lower than 
the free energy of the normal state (global criteria; 0.707 = l/\/2 is Chandrasekhar- 
Clogston limit of breaking BCS superfluidity at weak coupling). Right panel. Unstable 
regions according to criteria 2 and 3: superfluid solution should be a minimum of the 
free energy (local criteria), pressure of the superfluid state is higher than the pressure of 
normal state (global criteria). Requirement 2 excludes states above doted line. Criteria 
2 is the strongest among local criteria, and it leaves gapless states only with one Fermi 
surface. Local minimum of free energy (C > 0) is equivalent that the number susceptiblity 
matrix is positive definite. 



A, B, C, D, E involve Matsubara sums over bosonic frequencies 
ujn = (2n + 1)t:T and 3-d integration, which is done numerically. Ana- 
lytic expressions can be obtained in weak coupling, using BCS hierarchy 
of scales, <5/i, A << /.t (integration in thin shell around /i). Concluding, in 
the weak coupling gapless phases, (5/i > A, are unstable (for both Gold- 
stone and Higgs mode i? < 0, i? < 0, C < 0). In the weak coupling, 
the first order phase transition superfluid to normal state takes place, at 
(5/x = A/V2, before the gapless state developes, <5/i = A; i.e. for the window 
A/-\/2 < (5/i < A Higgs and Goldstone modes are fluctuations around the 
meta-stable solution. To find stable gapless phases one should explore the 
strong coupling regime. 



3. Analysis of stability at T = 0. 

Using the quasiparticle dispersion law, we identify regions in (J/i, /i) where 
the system has gapless excitations, i.e. momenta where energy is zero. Gap- 
less fermionic excitations live at one and two spherical surfaces in momen- 
tum space. We map regions with gapless excitations and regions where 
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Fig. 2. The four dot-dashed hnes show how fi varies as function of 5fi for four differ- 
ent values of the dimensionless variable k = 7r/2\/2mAas. Values k << —1, < 0, 
correspond to being deep in the BCS regime, while k >> 1, > 0, deep in the BEC 
regime. Region between curves 2 and 3, local and global minimum of free energy, corre- 
spond to meta-stable gapless state. Moving along k = 1.71 (BEC), point P correspond 
to the largest Sfj, where superfluidity is globally stable, S/i/A 1.59; after point Q, 
5fj,/A 1.66, it is locally unstable. Curves 2 and 3 run closer as 5fi increases, and con- 
verge asymptotically for Sfi >> A. Higgs mass is zero along curve 2. Hence Higgs mass 
is small along curve 3 and gets smaller as Sfi increases. 

coefficients A, B etc. are negative, indicating instabilities to the (6fi, ^) 
space. 

Stabifity criteria are 

1. Real speed of sound for Goldstone and Higgs modes {B, E > 0), 

2. SF is a local minimum (C > 0), 

3. SF is a global minimum (17^ — < 0), 

which are equivalent to 

1. Real speed of sound is equivalent to positive screening masses, 

2. Local minimum is equivalent to positive number susceptibility. 

At the Fig. 1, the ((5/i, fi) plane is divided into three regions correspond- 
ing to different number of gapless surfaces in momentum space; also unsta- 
ble regions according to criteria 1 and 3 are shown. 

Requirement 1 renders all states above the dotted and dashed lines 
unstable, which leaves gapless states with one and two Fermi surfaces as 
possiblity. Requirement 3 renders all states above solid line unstable. 

At the Fig. 1, unstable regions according to criteria 2 and 3 are shown. 
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Fig. 3. Higgs mass, = C/D, along k = 1.71 (g = l/kpa^ ~ 1.31) curve as 5fi 

increases. For small 5fi, the Higgs mass is constant in the gapped superfluid phase. 
Indeed, in the gapped phase free energy of superfluid is independent of 5fi, superfluidity 
is favored over normal and it is locaUy stable [C > 0), hence m^f is independent of 
5fi. In gapless phase, Higgs mass decreases. Superfluidity wins over the normal, hence 
the smallest mu in this regime (C describes oscillations about the global minimum). 
Between P and Q, there is meta-stable region; below Q, it is locally unstable (C < 0, 
mu is imaginary). From Sfj, = to P, mu drops by a factor of 7.5. 

However, we still learn from criteria 1 about the tendency of the sys- 
tem to go to another state. Real speed of sound for the Goldstone mode, 
which is equivalent that Meissner mass is real (or current instability), is 
more stringent than the real speed of sound of Higgs mode. Both criteria 
point towards a non-homogeneous phase. We learn, that LOFF like state 
is possible only in weak coupling at the BCS side. 

The global minimun of free energy is the strongest requirement. There 
is a sliver of parameter space corresponding to meta-stable superfluid state 
(it is a local minimum, according to 2, but not a global one of free energy, 
according to 3). Gapless superfluid state with one spherical Fermi surface 
is at BEG side (^ < 0). 



4. Parameters of the Higgs Lagrangian 

To be more concrete, we solve the gap equation for various scattering 
lengths and see where we land in the parameter space. At the Fig. 2, the 
four dot-dashed lines show how /i varies as function of (5/i for four differ- 
ent values of the dimensionless variable k — tt Region between 
curves 2 and 3, local and global minimum of free energy, correspond to 
meta-stable gapless state. Moving along k = 1.71 (BEG), point P corre- 
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Fig. 4. Numerical solution of the ODE, Eq. 8, for the condensate in a vortex configu- 
ration. Vortex center is at r = 0. A(r) increases from to A, as we go away from the 
inner core of the vortex, over a length scale rg = y/E/3C. Units are in A at both axis. 

spond to the largest Sfi where superfluidity is globally stable, (5/i/A « 1.59; 
after point Q, J/i/A « 1.66, it is locally unstable. Higgs mass is zero along 
curve 2. Hence Higgs mass is small along curve 3 and gets smaller as Sfj, 
increases. Thus, there is a light Higgs in the gapless BEC superfluid state. 

At the Fig. 3, we depict Higgs mass, mf^ — C/D, along n — 1.71 curve 
as (5/i increases. For small (5/i, the Higgs mass is constant in the gapped 
superfluid phase. In gapless phase, Higgs mass decreases. Between P and 
Q, there is meta-stable region; below Q, it is locally unstable (C < 0, mn 
is imaginary). From Sfi = to P, mn drops by a factor of 7.5. 

At point P: 

a. At P, there is the smallest Higgs mass in the regime where superfluidity 
is favored over normal state. 

b. Note about reliability of mean field. At point P, g = 1.31 and as k in- 
creases, g increases even further. In BCS g < (large), BEC g > (large), 
since g — l/kpUs large g corresponds to small Os, which means that the 
mean field is reliable. 

c. Effective theory near P: gapless fermions living on one spherical surface 
in momentum space, massless fluctuations in the phase of the condensate 
and massive but light fluctuations in magnitude of condensate. It is inter- 
esting to probe the spectrum. 

d. Quantum corrections may alter small Higgs mass (beyond mean field). 

It is interesting to explore part c, but even before a detailed study there 
is another striking consequence of our results. 
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Fig. 5. Square of vortex radius as a function of a distance R from the center of a trap 
(R ~ Sfi). Two slightly different curves correspond to two ways of calculation. Bottom, 
smooth curve, first 3-d momentum integrals are taken numerically followed by taking 
a limit of small momentum. Top, curve with a casp, the order of integration and small 
momentum limiting procedure is interchanged, which is apparantly wrong. 



5. Vortex in a trap 

The correlation length tq , or the typical length scale at which the magnitude 
of condensate A varies in field configurations that arise when the system 
is excited, is inversely proportional to the mass of the Higgs mode of the 
system, tq I/jtih- For example, rp governs the size of the outer core of 
a vortex configuration in a superfluid phase. This can be seen from the 
classical field equations for A(r) = f {r)e^'^'^'^\ f{r) = A + p{r), for static 
configuration, where A is the ground state value and p{r) is fluctuation of 
the condensate, (r, 0) are cylindrical polar coordinates, and the vortex is 
at r=0, 

A(r) - roV^A(r) = const 

/(r^0) = 0, /(r->oo) = A (8) 

ODE is solved numerically, Fig. 4, with boundary conditions: zero conden- 
sate at the center of the core, and it tends to a constant ground value at 
infinite r. For a vortex configuration, winds around by a multiple of 27r as 
we traverse a loop around the vortex. A(r) increases from to A, as we go 
away from the inner core of the vortex, over a length scale rg — yE/iC . 
The fact that C (Higgs mass) is numerically small close to the point P in 
parameter space, will manifest itself in an increased size for the outer core 
of the vortex. 
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To construct an actual vortex solution, it will be important to include 
the fourth order term in the effective action. However, the coefficient of this 
term is dimensionless and does not introduce a new scale, hence our basic 
arguments remain valid. 

To sec the effect quantitaivcly, we plot square of vortex radius as a 
function of a distance R from the center of a trap (harmonic trap in BEC; 
distance from the center of the trap is proportional to mismatch in chemial 
potentials, R ^ At the Fig. 5, in the gapped region the outer radius of 
the vortex is constant, and in the gapless region it increases monotonically 
on increasing R (or Sji). rg formally diverges at some point in the exterior 
of the trap, which corresponds to the transition from superfluid to normal 
state. In reality this divergency is absent, since there are no vortices in the 
normal state, i.e. our formula does not apply then. Two slightly different 
curves correspond to two ways of calculation. At the bottom smooth curve, 
3-d momentiim integrals arc taken numerically first followed by taking a 
limit of small momentum. At the top curve with a casp, the order of inte- 
gration and small momentum limiting procedure is interchanged, which is 
apparantly wrong. Here the long wavelength limit is not well defined, hence 
the order of operations matters. 

We estimate ro from the fact that it is the length scale at which the 
condensation energy is comparable to the kinetic energy of a superfluid el- 
ement close to the vortex, Ek — Econd- Taking into account the velocity of 
superfluid matter near a vortex, which diverges towards the center of a vor- 
tex V = 2^^6^ we estimate the kinetic energy, E^ = n^^. Condensation 
energy is proportional to a condensation energy density, Econd = necond- 
We get for the vortex radius ro ^ ^/^Scond- Therefore, in the gapped phase, 
ro is constant, since condensation energy density does not change; in the 
gapless phase condensation energy density decreases with increasing S/j,, 
hence ro increases steeply in the gapless phase; ro diverges at the transi- 
tion from superfluid to normal state, because this formula does not work 
in the normal phase (there are no vortices). In conclusion, ro continuously 
increases in the gapless phase. Vortex size grows steeply by entering into 
the gapless phase. 

We suggest to tune parameters of the trap, the number of particles of 
the two species, A''i and N2, and the scattering length, Og, so that there is a 
sufficiently wide region in position space where the atomic system is in the 
gapless BEC phase. We predict, that effect of sharp increase in vortex size 
can be seen. Probably, it requires very flat traps to realize this phenomena 
in a wide enough region to be observed. 
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6. Conclusion 

We obtained an effective theory for collective modes, Goldstone and Higgs, 
for the system of interacting two fermionic species. There is an interest- 
ing low energy content of effective theory, which includes gapless fermions, 
massless Goldstone mode, and very light Higgs mode. Prom our phase dia- 
gram, gapless state is stable at the BEG side of the resonance, with gapless 
fermions residing on one spherical Fermi surface in momentum space. Insta- 
bility towards non-homogeneous LOFF state occurs only at weak coupling, 
i.e. LOFF is favored away from the resonance. 

We would like to study the low energy effective theory further. Apart 
from that, we suggest a possible experimental evidance/signature of the 
gapless phases. We predict a sharp inc;reasc in outer core of a vortex, i.e. 
vortex size, upon entering into the gapless phase. To observe this dramatic 
effect will require tuning of the parameters of the trap to a gapless BEG 
phase. 

It will be interesting to study the core structure of the vortex, includ- 
ing inner and outer core, as well as interaction between two vortices. This 
requires solving Bogoluibov-De Ginnas equation and/or constructing the 
Landau- Ginsburg functional to the fourth order. However, we believe that 
our prediction with respect to vortex size will not change qualitatively. Ac- 
cording to our estimate, the mean field treatment is reliable in the region 
of interest. Probably, the main improvement will come from incorporating 
momentum/energy dependent Fermi interaction. 
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